Orbital-dependent two-band superconductivity in MgB2 
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We show that a two-band model with k-dependent superconducting gaps well describes the trans- 
mission and optical conductivity measured for MgB2 thin films. It is also shown that the two-band 
anisotropic model consistently describes the specific- heat jump and thermodynamic critical magnetic 
field He- A single-gap anisotropic model is shown to be insufficient to understand consistently optical 
and thermodynamic behaviors. In our model, the pairing symmetry in each band has an anisotropic char- 
acteristic which is determined almost uniquely; the superconducting gap in the cr-band has anisotropy 
in the ab-plane and the gap in the 7r-band has a prolate form exhibiting anisotropy in the c-direction. 
The anisotropy in the cr-band produces rather small effects on the physical properties compared to the 
anisotropy in the vr-band. 



Recently discovered superconductor MgB2 with a rel- 
atively high Tc=39K has attracted much attention in 
condensed-matter physics. ""^ An s-wave superconduc- 
tivity (SC) was soon established by experiments, e.g., 
a coherence peak in ^^B nuclear relaxation rate^ and 
its exponential dependence at low temperatures."^'^ An 
isotope effect has suggested phonon-mediated s-wave 
superconductivity.^ In contrast to its standard proper- 
ties, there have been several reports indicating unusual 
properties of the superconductivity of MgB2. Several 
studies reported two different superconducting gaps: a 
gap much smaller than the expected BCS value and that 
comparable to the BCS value given by 2A = S.SSfcsTc- 
Their ratio is estimated to be Amin/^max ~ 0.3 — 0.4 
using several experiments. It is also reported that 
the specific-heat jump and the critical magnetic field are 
reduced compared to the s-wave BCS theory.^'^° The 
other typical characteristic is a strongly anisotropic up- 
per critical field in c-axis-oriented MgB2 films and single 
crystals of MgB2.""^^ 

The unusual properties of MgB2 suggest an anisotropic 
s-wave superconductivity or a two-band superconductiv- 
ity. The band structure calculations predicted multi- 
bands originating from a(2px,y) and 7r(2pz) bands. In 
the ARPES measurements performed in single crystals of 
MgB2 three distinct dispersions approaching the Fermi 
energy were reported. 

There have been several studies on the anisotropy of a 
superconducting gap.^^"^° The two-gap model is shown 
to consistently describe the specific heat^^'^° and the up- 
per critical field Hc2^^ with the adoption of the effec- 
tive mass approach. In this paper, we examine opti- 
cal properties and thermodynamics to determine the k- 
dependence of the gaps. The k-dependence of two gaps 
in MgB2 has never been investigated thus far. We show 
that the optical transmittance, conductivity, specific-heat 
jump, and thermodynamic critical field are well de- 



scribed by a two-band superconductor model with differ- 
ent anisotropics in k-space. The symmetry in k-space 
is determined almost uniquely according to these exper- 
iments. First, we show that the single-gap model is in- 
sufficient to understand consistently optical and thermo- 
dynamic behaviors. Second, the two-gap model with dif- 
ferent symmetries in k-space is shown to be sufficient to 
understand optical and thermodynamic properties. 

The optical conductivity for anisotropic s-wave SC is 
investigated and compared with available data for MgB2. 
A simple angle-dependent generalization of the Mattis- 
Bardeen formula^^ is used to calculate the optical con- 



ductivity. The density of states N{e) 



is generalized to N{t) — (Ree/, /e^ — A? where the 



bracket indicates the average over the Fermi surface. We 
employed the following formula at T=0: 



CTls(w) 



= - / dE[N{E)N{uj 
^ Jo 

-(Re ^ ^^ =)k{Re 



E) 



A. 



E^ 



o'2s/ci„ is similarly generalized. The anisotropic order 
parameters considered in this paper are: 



Aci(k) = A(l + acos(2e')), 
A,2(k) = A(l-6cos2(0)), 
Aah(k) ^ A(l + ccos(6(/))). 



(2) 
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Here, 9 and (j) are the angles in the polar coordinate where 
is the polar angle with respect to the c-axis. The pa- 
rameters a, b and c determine the anisotropy. Ad is 
a prolate form gap for a > and is oblate for a < 0. 
^c2 {i> > 0) shows the same anisotropy as Ad for a < 0. 
Aab indicates a cylindrical gap with anisotropy in the a-b 
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FIG. 1. Transmission Ts/Tn for anisotropic s-wave models. 
The solid curve without marks is the Mattis-Bardeen result. 
Squares, triangles and circles are for the prolate, ab-plane 
anisotropic, and oblate gaps, respectively. The oblate form 
shows a small increase compared to other types. For the oblate 
gap, Ao = Amai, and for other types, Ao = A. 

plane corresponding to the a band. The integral in eq.(l) 
is detemined numerically by writing the average over the 
Fermi surface with elliptic functions. 

First, wc show that the one-band model is insufficient 
to understand consistently optical and thermodynamic 
behaviors. In Fig. I the transmission T5 at T = is 
shown as a function of the frequency lo. The following 
phenomenological expression for Ts/Tn is employed to 
determine the transmission curve theoretically,^^'^'* 

75 
Tn 

1 

[Tlf' + (1 - r^/^)(ai/a„)]2 + [(1 - Tl/')ia,/a^W ' 

(5) 

where cti and a2 are real and imaginary parts of the op- 
tical conductivity, respectively. T/v is determined from 
the expression for the ratio of the power transmitted 
with a film to that transmitted without a film given as 
Tn = 1/[1 + c7ndZo/{n + 1)]^. Here, d is the film thick- 
ness, 71 is the index of refraction of the substrate, and 
Zq is the impedance of free space. We have assigned the 
following values: d = I0~^ cm, n ^ 3, Zq — 377^1, and 
cr„ w 8 X lO^f^-icm-i. Then we obtain Tn ^ 0.014. The 
theoretical curves for Ts/Tn are shown in Fig.l; they 
have peaks near oj ~ 2Ao. For the oblate, its peak shows 
an increase only twice the normal state value, while the 
prolate and ab-plane anisotropics show more than twofold 
increases. The experiments show an approximately 2.5- 
fold increase^^ which supports the prolate or ab-plane 
anisotropic symmetry. However, the temperature depen- 
dence of the ratio H^2/H^2' which increases as the tem- 
perature decreases, ^'^ indicates that A(k) has an oblate 
form instead of a prolate form^^ in contrast to Ts/Tn- It 





0.5 1 1.5 2 2.5 3 
Q) X 4meV 

FIG. 2. Transmission Ts/Tn for the two-band model. The 
data points (large solid circles) are taken from ref.22 at 
T = 4K. The solid line without marks shows the Mat- 
tis-Bardeen result with 2A — 5meV. The dotted line is for the 
single gap of prolate type with 2Amax = 9meV and a — 0.5. 
Triangles are for the single-gap model of oblate form with 
2Amaa; = 9meV and a — 0.33. Others are for the two-band 
gap model where the ab-plane anisotropic form for the cr-band 
and the prolate form for the vr-band are assumed. The param- 
eters are the following. Solid circles: 2Amax = 8.5meV and 
c — 0.33 ((T-band: weight 0.45); 2Amax ~ 6meV and a = 0.33 
(TT-band: weight 0.55). Open circles: 2Amax ~ 9meV and 
c — 0.33 ((T-band: weight 0.45); 2Amax ~ 6meV and a = 0.33 
(TT-band: weight 0.55). Squares: 2Amax ~ lOmeV and c = 0.5 
((T-band: weight 0.4); 2Amax ~ 7.5meV and a — 0.5 (yr-band: 
weight 0.6). 

is also difhcult to describe the thermodynamic quantities 
such as the specific-heat jump at T = Tc and the thermo- 
dynamic critical magnetic field He within the single-gap 
model consistently. The specific-heat jump at Tc is given 
by 

AC(re) ^ J2_ (zy 

jcTc 7C(3) (z4) ' 

where 7c is the specific-heat coefficient and z is an 
anisotropic factor of the gap function, (z") is the average 
of z" over the Fermi surface. The experiments indicate 
that this value is in the range of 0.76~ 0.92;'^'-'^° the fit- 
ting parameters must be a ~ 0.3, —a 0.5 and c 0.3 
for the prolate, oblate and ab-anisotropic types, respec- 
tively. We must assign different values to parameters a 
and c in order to explain the thermodynamic critical mag- 
netic field i?c- The ratio of Hc{T = 0) to the BCS value 
is given as 

H,{T = Of/-fcT! - (67r/e2^)(z^) = 5.94(^2). (7) 

Thus to be consistent with the experimental results, 
(z^) should be less than 1; a should be small, a ~ 0.07, for 
the prolate form, and the ab-plane anisotropic and oblate 
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forms (a < 0) are ruled out since (z^) > 1. In Table I, 
wc siimmarize the status for the single-gap anisotropic 
s-wave model applied to MgB2. As shown here, it is 
difficult to understand the physical behaviors measured 
using several experimental methods consistently within 
the single-gap model. 

Here, a two-band model with two different anisotropies 
is investigated. We assume that the hybridization be- 
tween a and tt bands is negligible that the optical con- 
ductivity is given by 

Cr = ■WcrCTcr + W^ttCTtt, (8) 

where and (Ttt denote the contributions from a- and 
TT-bands, respectively. The transmission Tg/Tpf in Fig.2 
shows that the theoretical curve is in good agreement 
with the experimental curve. The effect of a&-anisotropy 
for the transmission Ts/Tm is small. The optical con- 
ductivity is also described well by the two-band model 
as shown in Fig. 3. We assign the following parameters 
to the best fit model in Figs. 2 and 3. The cr-band has 
ab-plane anisotropy with c « or less than 0.33 and the tt- 
band has the prolate form gap (cigar type) with a « 0.33. 
The ratio of the weight of the cr-band to that of the tt- 
band is 0.45/0.55, which agrees with penetration depth"^ 
and band structure calculations.^^ The ratio of the min- 
imum gap to the and maximum gap is 0.35, which is in 

1.2 
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FIG. 4. Thermodynamic critical magnetic field Hc{T) nor- 
malized by Hc{T = 0)bcs- The bold dashed curve indicates 
data from rcf.lO and the bold solid curve indicates those ob- 
tained using the present two-band anisotropic model. The 
thin solid curve indicates the BCS results. The results for the 
cr-band (ab-plane anisotropic) and the 7r-band (prolate form) 
are also shown. 
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FIG. 3. Real part of the optical conductivity for the two-band 
anisotropic model. The data points (large solid circles) 
are taken from ref.22 The parameters for solid circles are 
'iAmax = 8.5meV and c = 0.33 (cr-band: weight 0.45); 
2Amax = 6meV and a = 0.33 (7r-band: weight 0.55). The 
parameters for open circles are 2Amax ~ lOmcV and c = 0.5 
(cr-band: weight 0.4); 2Amax = 7.5meV and a — 0.5 (vr-band: 
weight 0.6). The dashed line indicates the results obtained 
using the Mattis-Bardeen formula with 2A = 5meV. The 
dash-dotted line denotes the conductivity for the d-wave 



the range of previously reported experimental values.®'^ 
Let us mention here that the two-band isotropic model 
(z — I) describes inconsistently the observed (Tis/o"i„ 
since a shoulder-like structure is predicted in the two-gap 
isotropic model if the two gaps have different magnitudes. 
In Fig. 4 the thermodynamic critical magnetic field Hc{T) 
is shown for the single-band and two-band models with 
available data.^° We have simply assumed that the total 
free energy is given by the sum of two contributions from 
a- and 7r-bands: O = w^i^cr + Wj^il^. The experimental 
behavior is well explained by the two-band anisotropic 
model using the same parameters as those for Ts/Tn 
and (Tis/ain- We show several characteristic values ob- 
tained from the two-band model in Table II. Results of 
analyses of Hc2 and specific heat using the effective mass 
approach are consistent with those obtained using the 
two-band model. It has been reported that the 
increasing nature of H^2 1^0.2 with decreasing tempera- 
ture is explained by the two-Fermi surface model. The 
spccific-hciat coefhcient 7 in magnetic fields seems consis- 
tent with that of the multiband superconductor. ^^'^"^ 

We have examined the transmittance, optical con- 
ductivity, specific-heat jump and thermodynamic criti- 
cal magnetic field He of MgB2 based on the two-band 
anisotropic s-wave pairing model. We have shown that 
the two-gap model with different symmetries in k-space 
can explain the experimental results consistently. 
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TABLE I. Anisotropic parameters in the SC gap function used to fit several physical quantities. The upper four rows are for 

the single-SC gap model and the last row is for the two-band anisotropic model for comparison. The cross indicates that we 
cannot fit experimental data by the corresponding z factor. A in the column Hc2 indicates that experiments are explained 
qualitatively but not quantitatively.^^ The big circle show that we can fit the data using the same parameters in the column 
cTi. For the two-gap model the anisotropy parameter in the a-band can be small since the cr-band anisotropy produces only 
small effects on physical quantities. 
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TABLE IL Several physical quantities obtained by the two-band model with c ~ 0.33 (cr band) and o ~ 0.33 (tt band). 





"V/"'7r 


A„„„/A,„„,, 


Arc;:.) 


11, (n) 


Two-band 


0.45/0.55 


- 0.35 


- 0.82 


- 0.95 


Exp. 


0.45/0.55 


0.3-0.4 


0.76 - 0.92 
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